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Abstract 

These notes contain a rapid overview of the methods and results obtained in the field of 
propagation of waves in disordered media. The case of Schrodinger and Helmholtz equations 
are considered that describe respectively electrons in metals and scalar electromagnetic waves. 
The assumptions on the nature of disorder are discussed and perturbation methods in the weak 
disorder limit are presented. A central quantity, namely the probability of quantum diffusion 
is defined and calculated in the same limit. It is then shown that several relevant physical 
quantities are related to the return probability. Examples are provided to substantiate this, 
which include the average electrical conductivity, its fluctuations, the average albedo and 
spectral correlations. 

1 Introduction 

The study of wave propagation in random media gave rise to a huge amount of work especially dur- 
ing the last twenty years. Today, this field is split in two main subfields. One is concerned with the 
interplay between coherence and disorder in metallic systems and the second deals with the same 
problematics but for electromagnetic waves (this includes as well sound or gravity waves). Each 
one of the subfields has its own specificities and advantages, such as the effect of magnetic fields 
on transport and thermodynamics, interactions for metals and angular structure (spectroscopy) 
for waves. Unfortunately, the split between these two subfields has been growing so that their 
interplay is now certainly too weak in spite of the existence of a number of excellent reviews. 
These notes represent a very preliminary step towards a unified presentation of this field jjj . We 
have tried to present a general formalism that can apply to both situations (electrons and waves). 
It is centered around the existence of a basic quantity the probability of quantum diffusion which 
allows to describe either weak localization effects of the electrical conductivity of metals, spectral 
quantities of isolated electronic systems or the coherent and incoherent albedo, dynamical effects 
in multiple scattering of light by suspensions etc. 

The next three sections contain basic definitions and generalities on the kind of waves we 
consider and the model of disorder. The section 5 deals with the definition and the calculation of 
the probability of quantum diffusion in multiple scattering using various approximations. Then, 
the subsequent sections apply these results to a selection of examples taken either from electronic 
systems (conductivity and spectral correlations) or multiple scattering of electromagnetic waves 
(average albedo in optical systems). 



2 Models for the disordered potential 

We shall consider mainly two problems where waves and disorder are involved. The first one 
corresponds to the study of spinless electrons in disordered metals or semiconductors. For a 
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degenerate gas of free electrons of mass m and charge — e, the Schrodinger equation is governed 
by the Hamiltonian 

h 2 ip 

= ~2^ (V + ft A ) + V{Y)lp{Y) ■ (1) 

where B = V x A is the magnetic field. Effects associated to the band structure or to interactions 
between the electrons in the framework of Fermi liquid theory are accounted for by the replace- 
ment of the mass m by an effective mass. The potential V(r) decribes both the scattering by 
inhomogeneities and the confinement potential. 

The second problem we shall consider is the propagation of electromagnetic waves of frequency 
uo in the scalar approximation. The behaviour of the electric field ip(r) is obtained from the 
Hclmholtz equation 

- A^(r) - fc 2 /x(rMr) = k&(r) (2) 

where the function /z(r) = Se/e is the relative fluctuation of the dielectric constant, kg = and 
n is the average optical index. Under this form, the Hclmholtz equation has a structure similar to 
the Schrodinger equation and the waves are scattered by the fluctuations of the dielectric function. 
It is nevertheless interesting to notice that in the latter case, the strength \i of the potential is 
multiplied by the frequency ui 2 so that in contrast with electronic systems, a decrease of the 
frequency u> leads to a weaker effect of the disorder. 

To account for the effects of the disorder in cither case, we shall consider a random continuous 
function V(r), of zero spatial average, (V(r)) = 0, where (• • •) represents the spatial disorder 
average. The potential is characterized by its correlations especially by the two-point correlation 
function 

(V(r)V(r>))=B(r-r>) (3) 

For a potential V(r) which is localized enough, there exits a length r c which describes the fall 
off of B(r - r'). In the limit where the wavelength A ^> r c , the scattering events are statistically 
independent and we can consider the limiting case 

(V(r)V(r'))= 7 6(r-r') (4) 

Such a potential is usually called a white noise. This is the case we shall consider throughout 
these notes. For the case of the Helmholtz equation, the potential is taken to be V(r) = k 2 ^i(r) 
so that 7, in that case, has the dimensions of the inverse of a length. 



3 Perturbation theory for multiple scattering 
3.1 Single scattering - elastic scattering time 

For a very dilute system, we can first assume that a given incident wave of wavevector k is scattered 
only once into a state k' before leaving the system. The lifetime Tk of the state k is then given 
(using the notations of quantum mechanics) by the Fermi golden rule 



For the white noise case, we obtain 



— = -rPol (6) 

where p is the density of states per unit volume. To this time, we can associate a length, the 
clastic mean free path l e , defined as l e — vr e by using the group velocity v. But the Fermi golden 
rule is valid for short times t <§C r e and in perturbation with the scattering potential V. In order 
to go beyond these limitations, we need to resort to a more powerful tool namely the formalism 
of the Green functions. 
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The Green function associated to a wave equation can be defined as the response to a pulse i.e. 
to a ^-function perturbation. For the case of the hamiltonian (Q), the Green function G R,A (ri, r, e) 
is solution of 

(e-H± iO) G R ' A (n,T, e) = 5(r - r<) (7) 

We can define as well the free Green function Go(r,,r, e) in the absence of scattering potential, 
which is the solution of 

e + ^A r ± i0\ G*' A (r l: r, e) = S(r - n) (8) 
so that the Green function G R ' A {vi, r, e) can be also expressed as a solution of the integral equation 
G(r u r, e) = G (r 4 , r, e) + J G(r„ r', e)V(r')G (r' , r, e)rfr' (9) 
the solution of (g) is given (for d=3) by 

\ m e ±lkR 

G ' (Ti,r,e) = - 2 (10) 

27T?l ii 

with e = h 2 k 2 /2m. 

For the Helmholtz equation (||), we obtain a sequence of similar equations, but attention needs 
to be paid to the fact that the dispersion of the waves is now linear instead of quadratic for the 
Schrodinger case. The Green equation is 

(A r + fc 2 (l + M (r))) G(r h r, k ) = 5(r - r') (11) 

while the solution of the free Green equation obtained for /x(r) = 0, is 

i „±ik R 

G R ' A (r,,r,fco) = --— — (12) 
3.2 Electromagnetic waves 

The formalism of Green functions provides an appropriate technical framework for the study of 
solutions of wave equations for free systems namely without sources. It becomes essential for the 
study of the propagation of electromagnetic waves from a distribution of sources j(r) of the field. 
This is the problem of radiative transfer (2). For a pointlike source, we are back to the previous 
problem of the Green equation. In the general case, The Helmholtz equation (|J) needs to be 
replaced by 

AVXr) + fc 2 (l + /i(r))^(r) - j(r) (13) 

where ?p(r) is indeed a Green function i.e. it depends on the distribution of sources j(r). The 
equation (Oh can also be written in the form of the integral equation 



lfr(r) = / dr 4 j(r l )G (r 4 , r, fc ) - K J drV(r>(r')G (r', r, fc ) (14) 
which allows to consider separately the effects of the source and of the random potential fi(r). 

4 Multiple scattering expansion 

Either the expressions (||) or (|l4|) provide the starting point for a systematic expansion of the 
Green function in terms of the free Green function. It can be written 

G(r,r') = Go(r,r')+ / driGo(r,ri)y(n)G (ri,r / ) 

+ y"dr 1 dr 2 Go(r,r 1 )F(r 1 )Go(r 1 ,r 2 )^(r 2 )G (r 2 ,r') + ... (15) 
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We can now calculate the average Green function using the white noise potential (Q) . All the odd 
terms in the potential V disappear from (JL5|), and it remains 



G(r,r') =G (r,r')+7 / driGo(r, n)G (ri, ri)G (ri, r') + . . . (16) 

where we denote from now on the disorder average by 777 . By averaging over the disorder, the 
medium becomes again translational invariant and the Green function G(r,r') = G(r — r'). The 
average over the disorder generates all possible diagrams. Among them, there is a subclass called 
irreducible diagrams that cannot be split into two already existing diagrams without cutting an 
impurity line. It is possible to rewrite the average Green function or its Fourier transform in terms 
of the contribution of these diagrams only. We then obtain the so called Dyson equation 

G(k) = G (k) + G (k)£(k, e)G(k) (17) 

where the function £(k, e) is called the self-energy. It should be emphasized that, although the self- 
energy contains only the irreducible diagrams, there is an infinity of them. Thus, the calculation 
of £ is a difficult problem. For the white noise potential (Q), £ can be expanded in powers of 7. 
To first order, for the Schrodinger equation, we obtain 

^(vj^E^w ( 18 ) 

q 

where f2 is the volume of the system. The real part of £ corresponds to an irrelevant shift of the 
origin of the energies that we shall ignore. The imaginary part is 

ImEf(k,c) = -7rpo(e)7 (19) 
while for the Helmholtz equation, it is 

Im£f(k,e)=-^ (20) 

We emphasize again that the difference between these two expressions results from the two distinct 
dispersions of respectively the Schrodinger and Helmholtz equations. Higher orders terms in the 
expansion of £ are proportional to Im£{* times some power of the dimensionless parameter t4-. 
The contribution £1 describes the multiple scattering of the wave as a series of independent 
effective collisions. The higher corrections include interference effects between those successive 
scattering events. The weak disorder limit kl e ^S> 1 amounts to neglecting these interferences. We 
deal then with the so called self-consistent Born approximation. We shall, from now on, consider 
only this limit. 

It is then a straightforward calculation to get an expression for the average Green function at 
this approximation: 

G RA (r u r, ko) = G^ A (r u r, k Q ) e^^ 21 ' (21) 

To conclude this section, we would like to notice that although this expression has been obtained 
for the case of an infinite system, this restrictive assumption can be released and we need to 
consider, for this relation to be valid, only systems of sizes L 3> l e - 



5 Probability of quantum diffusion 

The quantities of physical interest are usually not related to the average Green function but instead 
to the so called probability of quantum diffusion which describes the probability for a quantum 
particle (or a wave) to go from the point r to the point r' in a time t. Once we average over the 
disorder, we shall see that this probability P(r,r',t) contains mainly three contributions: 

• i. The probability to go from r to r' without scattering. 
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• ii. The probability to go from r to r' by an incoherent sequence of multiple scattering, which 
is called the diffuson. 

• iii. The probability to go from r to r' by a coherent multiple scattering sequence. We shall 
calculate one such coherent process, called the cooperon. 

We shall first define the probability for the Schrodinger case. Please notice that througout this 
section we shall take h — 1. To that purpose, we consider a gaussian wavepacket of energy eo. 
We shall also assume that around eo the density of states is constant. Then, we can write for the 
Fourier transform P(r, t 1 ,lS) of the probability P(r, r',i) the expression: 

1 



P(r, r', u) = G*(r, r', e )G A (r>, r, e - u) (22) 

This probability is normalized to unity which means that cither 

J P(r,r',t)dr' = 1 (23) 



or 



P(r,r',w)dr' = - (24) 

UJ 

5.1 Free propagation 

In the absence of disorder, the Green functions in j2^ ) take their free expression fljo| ) and it is 
straightforward to obtain for the three dimensional case 

where i? = |r' — r| and v being the group velocity. This probability is indeed normalized. 

5.2 Drude-Boltzmann approximation 

In the presence of disorder, we need, in order to calculate the probability, to evaluate the average 
of the product of the two Green functions that appear in (p2|). The simplest approximation is 
to replace the average by the product of the two averaged Green functions. Here again, since we 
have calculated in the weak disorder limit the expression of G, we obtain 

iu>R/v-R/l e 

so that 

/P (r,r»dr' = — %— (27) 

J 1 - lUT e 

At this approximation, the probability is not normalized, but instead 

[ P {r,v',t)dr' = e-^ (28) 

It is then clear that some part is missing in the probability. The Drude-Boltzmann approximation 
overlooks a large part of the probability; since after a time f, it predicts that the wavepacket 
disappears. 
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5.3 The diffuson 



There is another contribution associated to the multiple scattering which can be calculated in the 
weak disorder limit kl e 3> 1 in a semiclassical way. Using the description we obtained previously 
for the calculation of the average Green function, we can associate || to each possible sequence 
C, of independent effective collisions a complex amplitude A(r, r',C). Then, using a generalization 
of the Feynman path integral description, we can in principle write the Green function as a sum 
of such complex amplitudes. 

Then, in order to evaluate the product of two Green functions, we notice the following two 
points. 

• i. Due to the short range of the scattering potential, the set of scatterers entering in the 
sequences for both G R and G A must be identical. 

• ii. For the effective collisions, the mean distance between them is set by the elastic mean free 
path l e A. Therefore, if any two scattering sequences differ by even one collision event, 
the phase difference between the two complex amplitudes, which measures the difference of 
path lengths in units of A will be very large and then the corresponding probability will 
vanish on average. 

We shall therefore retain only contributions for which the corresponding probability P<i(r, r', u>) is 

Pz(r,r» = 2^/Gf(r,r 1 )Gf_ (1) (r 1 ,r)Gf(r 2) r')Gt tl ,(r' ) r 2 )x 

x T L ,{r ll r 2 )dr 1 dr 2 (29) 
It is made of two multiplicative contributions. The first one is 

G? (r, ri )Gf (r 2 , r')Gf_ Jn, r)Gf_Jr', r 2 ) . 

It describes the mean propagation between whatever two points r and r' in the medium and 
the first (ri) (respectively the last (r 2 )) collision event of the sequence of scattering events. The 
second contribution defines the quantity ^(r^ra) which we shall call the structure factor of the 
scattering medium. In a sense, it generalizes to the multiple scattering situation the usual two- 
point correlation function in the single scattering case. We now use once again the assumption of 
independent collisions in order to write for T u (ri, r 2 ) the integral equation 

r w (r 1 ,r 2 )=7*(ri-r 2 )+7 / Sf (r x , r)Gf_Jr, r x )r w (r, r 2 )dr (30) 

This equation can be solved exactly in some geometries. For the infinite three dimensional 
space, we can make use of the translational invariance and get for the structure factor the expres- 
sion 

r«(q) = z p 7 w (31) 

1 - P (q, W)/T e 

where Po(q, w) is the Fourier transform of ( p6|) and is given by arctan 1 ^ l i ^ JT with q = |q|. 
Then, the probability rewrites 

D / s o / \ P)(q, w)/r e 

P rf (q,w = P (q,w)- — ry- 32 

1 - P (q, UJ)/Te 

Using this expression of Pz, the normalization of the total probability P = Pq + P<j can be readily 
checked namely P(q = 0, u>) — ~. For the semi- infinite space with a point source, it is also possible 
to obtain a closed analytical expression for the probability using the Wiener-Hopf method. But 
beyond these two cases, for simple finite geometries, it is not possible to obtain the solution of 
(|29j) without resorting to numerical calculations. We are then led to look for some approximate 
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solutions. An excellent one is the diffusion approximation obtained for large times ( > t c and 
large spatial variations r>I e . It is obtained by expanding the structure factor under the form 

ri).V ri ] 2 r w (33) 



ra) (34) 

where the diffusion coefficient is D = 4— = ^v 2 T e . At this approximation, we have between P& 
and r w the following relation 

P rf (r,r»~2 W 2 r w (r,r') (35) 

so that Pd, as well, obeys a diffusion equation. It is interesting to check the validity of the diffusion 
approximation. For an infinite system, and for j- = 1 the relative correction between the exact 
solution and diffusion approximation is 0.085 while for j- — 2.5, it is less than 5.10~ 3 . 

5.4 The cooperon 

With the normalized expression of the probability we just have obtained, it seems that we fulfilled 
the demand of evaluating all the relevant processes that contribute to the probability. But it could 
be, and it is certainly the case, that there are many other contributions that sum up to zero. 

For instance, we may consider the possibility which corresponds to the product of Green 
functions such as we considered before, but where now the two identical trajectories are time 
reversed one from the other. It is clear that if these trajectories are closed on themselves, there is 
no phase difference left between them. This requires that the system has time-reversal invariance 
namely that G R ' A (r, r', t) — G RA (r' ', r, t). This relation does not hold anymore in the presence of 
a magnetic field for electronic systems. 

The contribution, we shall call P c , of this process to the total probability can be evaluated as 
we did before for the diffuson. Thus we have instead of (|29|) 

^c(r,r» = ^-|Gf(r,r 1 )Gf(r 2 ,r')Gf_ a; (r',r 1 )Gt w (r 2 ,r)x 

x r' w (r 1 ,T 2 )dr 1 dr 2 (36) 
where the new structure factor T'^ is solution of the integral equation 

r^n.ra) =7<y(n-r a )+ 7 J Gf(r i y)Gt^i,r'')K(r'',r- 2 )dr'' (37) 

ft j± 

Notice that unlike the structure factor of the diffuson, the new combination G e (ri, r")G e _ a ,(ri, r") 
cannot be simply written in terms of the probability Pq. But as before, we can evaluate P c (j, r', lu) 
in the diffusion approximation (i.e. for slow variations) and we obtain 

P c (r, r', u) ~ E^£l£l [ j G? (r, n)cf (r', n)dn] 2 (38) 

Since in the presence of time reversal invariance, we have 

Sffa.iOSffa.r) = Gf(r 1) r)Gf(r,r 1 ) (39) 
then, r^(ri,r 2 ) = r w (ri,r 2 ) and finally, 

P c (r, r', w) = P d (r, r, u)g(i - r') 2 (40) 



r w (r,r 2 ) = r w (ri,r a ) + (r - ri).V ri r w + - [(r - 



which together with the integral equation (|30|) gives 

[-iw-DA ri ]r w (ri,r a ) = ^^^i 

2-KpQTi 
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where in 3d we have the relation 

m - sj ^ e ~ R/2l ° (4i) 

and R = |r — r'|. For R = 0, i.e. for r = r', we have 

P c (r,r,w) =P d (v,v,Lj) (42) 

namely, the probability to come back to the initial point is twice the value given by the diffuson. 
The contribution of the cooperon P c to the total probability is given by 

[ P c (r,r',u)dr' = P d (r,r,u)^- (43) 

for any space dimensionality. How does this contribution compare with the diffuson contribution ? 
We have found that Pd{q = 0, t) ~ 1 while P c (q = 0, t) 2t ^ (Dt) 3 / 2 ^ or sman enough times. Then, 
P c (g = 0,i) is maximum for t ~ r e and given by P c (g = 0, r e ) = ^ } )d _ 1 . Thus, the contribution 
of P c to the total probability is vanishingly small for kl e ^> 1 and for a space dimensionality d > 2. 
But although it is very small, P c must be compensated by another contribution in order to restore 
the normalization of the probability. The additional contributions result from other irreducible 
diagrams. But the subsequent terms in this series are not known. 

We would like to conclude this section on the cooperon by emphasizing that although obeys 
a diffusion equation, it would be meaningless and incorrect to state that P c (r, r',w) obeys it as 
well. The exact statement is that for r = r', P c and P d are proportional and that P d obeys a 
diffusion equation. 



6 Radiative transfer 

6.1 Local intensity and correlation function 

In the previous sections, we have defined the quantum probability for electronic systems. The 
probability P is directly related to quantities that are physically measurable like the electrical 
conductivity or the magnetic response e.g. the magnetization ||. For the study of the prop- 
agation of electromagnetic waves in disordered media, the quantity which is usually measured is 
the local intensity of the field or its correlation function ||. As we discussed previousl y, a nd by 
definition of the Green function, the radiative solution (r) of the Helmholtz equation (|13|) with 
a localized source at point R = is "0c ( r ) — G e (0, r). The correlation function of the field is then 



^(r)V*_Jr') = Gf (0,r)G^(r',0) (44) 

It is not directly related to the probability of quantum diffusion P(r,r',u>). But the radiated 
intensity J(r) defined by 

4-7T 

7(r) = -\Mr)\ 2 

4-7T 

= — Gf(0,r)Gf(r,0) (45) 

is indeed related on average to the probability P and 

4tt 



7(r) = — G«(0,r)G^(r,0) (46) 
c 

Using for the probability the following relation which is the counterpart, for the Helmholtz equa- 
tion, of the relation ( |22] ) 

P d (v, r') = — G«(r,r')G^(r',r) (47) 
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we obtain J(r) = Pd(0, r). From now on, we shall denote I(r) the average intensity. 

We can rephrase what we did before in order to calculate the various contributions to the 
intensity that come respectively from the Drude-Boltzman, the diffuson and the cooperon approx- 
imations. The first contribution is given by 

= m 

It corresponds to the contribution to the radiative intensity of waves that did not experience any 
collision on a distance R from the source. 
The diffuson contribution is given by 

I d (r) = — J dr x dr 2 1 ^ e (r ! ) 1 2 r w=0 (r t , r 2 ) | Gf (r 2 , r ) | 2 (49) 

and finally, the contribution of the cooperon to the intensity is 

Ajr f j. d A 

7 c (r) = — / dr 1 dr 2 ^(r 1 ) 1 /. e (r 2 )r(r 1 ,r 2 )G e '(r 2 ,r)G e (r,n) (50) 
where we used the notation r w= o = r. In the diffuson approximation, the intensity Id rewrites 

7 d (r) = P d (0,r) = J-r(0,r) (51) 

and like Pd, it obeys the diffuson equation 

- DAI d (r) = S(r) (52) 

whose solution in the 3d free space is 

We shall now apply all the considerations developed in this section to the calculations of 
physical quantities in some specific situations both for metallic systems and for the propagation 
of electromagnetic waves in suspensions. 



7 Example 1. The electrical conductivity of a weakly dis- 
ordered metal 

We previously defined a weakly disordered metal as a non interacting and degenerate (spinless) 
electron gas (at T = 0), moving in the field of defects and impurities described by the white noise 
potential ([§). 

The average electrical conductivity <j(u>) calculated in the framework of the linear response 
theory JjJ is given by the Kubo formula: 

CTH = 2^Sfi £ k * k '* Gf(k,k')G^(k',k) (54) 

k,k' 

Then, we see from this definition, that the structure of the conductivity is up to the product 
k x k' x very similar to those of the quantum probability P. Therefore, and just as we did before, 
the very definition of cr(w) leads us to the following set of approximations. 
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7.1 The Drude-Boltzman approximation 

It is obtained by replacing the average of the product of two Green functions by the product of 
the averages, namely 

G«(k,k')G^(k',k) ~ cf (k, k')Gf_ w (k', k) (55) 
where the Fourier transform of the averaged Green functions (|2l]) is 

Gf'Vk') = Gf' A (k)5 k , k - = . R (56) 

e e(K) ± 

Then, we obtain for the conductivity cto(w) at this approximation the following expression 

a (u;) = — /"dr'P (r,r» (57) 

where Po is the quantum probability calculated at the same approximation (p6|). Using the ex- 
pression (p7|), we obtain 

— 2 

ff (w) = I^—^ (58) 

ml — iu>T e 

which is the well-known Drude expression. It must be noticed that because of the Kronecker delta 
function that appears in the average Green functions, the scalar product k x k' x reduces simply to 

k\ and eventually, after averaging, to -f where kp is the Fermi wavevector. 

7.2 The contributions of the diffuson and the cooperon 

Here again, we approximate the average product in the relation ([54]), using the same scheme we 
used for the diffuson. The contribution <Jd(w) of the diffuson to the conductivity is thus 

°aM = ./^o 1 a P d(0,")J2 k x k' x P (k,0,uj)P (k',0, u ) (59) 



k,k' 



where the function Po(k, q, to) = 2^cf (k + ^)G c _ w (k — |). Since the function -P (k, q = 0, uj) 
depends only on the modulus of the wavevector k and not on his direction, the angular integral in 
the previous expression gives a vanishing contribution namely cr^(iu) = 0. Then, it is interesting 
to notice that although the diffuson gives the main contribution to the quantum probability, its 
contribution to the conductivity which indeed measures such a probability, vanishes identically. 
We evaluate, the contribution a c (oj) of the cooperon using the relation ( |36| ) so that 

k Q 

the sum in the brackets is straightforward so that in the diffusion approximation, we obtain 

tr c (u) = -^-P 4 {T,T,u) (61) 
Using now the relation (^) between the cooperon and the diffuson, o~ c (uj) rewrites 

a e (u) = -— f dv'P c {T,v',uj) (62) 
m J 

and the total conductivity at this order is now given by 

a(u) = ^ / dr'(P (r,r',o;)-P c (r,r>)) (63) 

The conductivity <j(u>) is reduced by the coherent (cooperon) contribution. This correction is 
called the weak localization correction to the conductivity |^|, ^, |10"| . Relatively, this contribution 
of P c is much larger than the normalization correction to the total probability. This is because P c 
is now compared to Pq and not to Pd which represents the main contribution. 
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7.3 The recurrence time 



It is of some interest to study the dc conductivity <j(lo = 0) using a slightly different point of view 
HI. ^From the relation (j63|), and using the dc expression oo = " e m Te , we obtain for the relative 
correction to the conductivity, 

8(7 1 f°° 

— = =— / dtP c (r 7 r 7 t) (64) 

where po is the density of states per unit volume. 
Consider now the quantity Z(t) defined by 

Z(t) = [ P d (r,r,t)dr (65) 
Ja 

where the integral is over the volume CI of the system. It represents the return probability to 
a point r averaged over all those points. This quantity which characterizes the solutions of the 
diffusion equation is sometimes called the heat kernel in the literature. The time integral of Z(t) 
defines the characteristic time Tr 



T R = dr P d (r,r,t)dt = Z(t) dt (66) 
Jn Jo Jo 

called the recurrence time. It measures the space average of the time spent by a diffusive particle 
within each infinitesimal volume. Tr diverges, as stated by the Polya theorem ]TH , for a random 
walk in the free space of dimensionality d < 2. Then, we define a regularized expression for Tr 
given by the Laplace transform 

/>OC 

T R (s) = / Z(t)e~ st dt (67) 
Jo 

This expression of the recurrence time may be interpreted by saying that Tr(s) selects the con- 
tribution of all the diffusive trajectories of lengths smaller than L s = y/D/s. 

Then, the expression (|64| ) of the relative contribution of the cooperon to the conductivity 
rewrites 

c a /-oo A 

— = - — dtZ(t)e- st = —rT R (s) (68) 
cr ttTi J Q nh 

where the energy defined by A = j s the mean level spacing measured at the Fermi level 
between the energy levels of an electron gas confined in a box of volume CI. The weak localization 
correction to the dc conductivity can be essentially expressed in terms of a purely classical quan tity , 
namely the recurrence time Tr. When it diverges, we need to use its regularized version (|67|), 
where now the length L s can be given a physical meaning. The cooperon correction results from 
the absence of any relative phase between two time reversed multiple scattering trajectories. Any 
interaction of an electron with an external perturbation may destroy this phase coherence and 
then the contribution of the cooperon. For a large class of such perturbations which includes 
inelastic collisions at finite temperature, Coulomb interactions or coupling to other excitations, we 
can define a phenomenological length called the phase coherence length and a phase coherence 
time such that L 2 ^ — Dt^. Thus, is the length L s used to regularize the recurrence time. 
We shall see later other examples of physical quantities that can be expressed using either the 
heat kernel or the recurrence time. 



7.4 Conductance fluctuations- Spectral determinant 

We could as well calculate the so-called conductance fluctuations using the properties of the heat 
kernel. The conductance G is related to the conductivity a by the Ohm's law G = aL d ~ 2 . We 
define the dimensionless conductance g = ^G. Its fluctuation defined by (Sg 2 ) = {g 2 ) — (g) 2 can 
be expressed in terms of Z(t) through 

i o f°° 

(Sg 2 ) = — dttZ{t)e- st (69) 

T D Jo 
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where td = L 2 / D is the diffusion time. For a one dimensional system (in the sense of the diffusion 
equation), we have 

(Sg 2 ) = ^ (70) 

It is interesting to recover the last two results on the average conductivity (or conductance) and 
its fluctuation using a systematic expansion of the recurrence time. To that purpose, we define 
the spectral determinant S(s) by S(s) = det(—A + s) where A is the Laplacian operator. Then, 
by definition of Z(t), we have the relation 

T R (s) = / Z(t)e- St dt = — lnS(s) (71) 
Jo ds 

up to a regularization independent of the Laplace variable s. This relation is valid for all space 
dimensionality. Consider now as a working example |^| the case of a Id diffusive wire of length 
L. In order to describe the case of a perfect coupling of the wire to the reservoir, we demand 
Dirichlet boundary conditions for the diffusion equation. This is in contrast to the Schrodinger 
equation for which this choice corresponds to Neumann boundary conditions. Then, the spectral 
determinant S(s) of the diffusion equation can be readily calculated and it is given by 

S(x) = -^L (72) 
sinn y/x 



with x — std- The average conductance Sg deduced from ( p8|) and the fluctuation (Sg ) can be 
written generally as 

Sg = 2^-lnS(x)=~ (73) 

(Sg 2 ) = 12^1n5(x) = A ( 74 ) 

In the limit x — » 0, we can expand the spectral determinant namely hiS(x) ~ x ^a — § + and 
recover the previous result. 



8 Example 2. Multiple scattering of light: the albedo 

We shall now give an example of the use of the quantum probability taken from the multiple 
scattering of electromagnetic waves in disordered suspensions. To that purpose, consider the 
scattering medium as being the half-space z > 0. The other half-space is a free medium which 
contains both the sources of the waves and the detectors. We also assume that both incident 
and emergent waves are plane waves with respective wavevectors kj = fcij and k e = fcs e , where 
§i and s e are unit vectors. The waves experience only elastic scattering in the medium so that 
after a collision, only the direction s changes while the amplitude k = — remains constant. 
The reflection coefficient a(§i,s e ) for this geometry is called the albedo and is proportional to 
the intensity /(R, §i,s e ) emerging from the medium per unit surface and per unit solid angle 
measured at a point R at infinity. If Fi nc defines the flux of the incident beam (related to the 
incident intensity Ii nc by Fi nc = cSIi nc ) where S is the illuminated surface in the plane z = 0, 
then the albedo is given by 

R 2 c 

a(si,s e ) = — — /(R, §i, s e ) (75) 

r inc 

In order to calculate the intensity /(R, §i,s e ), we need to calculate first the radiative solutions of 
the Helmholtz equation at the Fraunhoffer approximation. It is given by 

^ (si,Se)= / 'drdrV fc ( §ir - § * r ')G(r,r> ) (76) 
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and the intensity is then I(R, §i, s e ) = |-0 Wo (si, s e )| 2 . Before averaging over the disorder, the albedo 
looks like a random pattern of bright and dark spots. This is a speckle pattern [D| whose 
long range correlations are a specific feature of multiple scattering. By averaging over the disorder 
(either liquid suspension [[l4| [is], |l6) or rotating solid sample |l7|), the speckle is washed out and 
the two surviving contributions are respectively associated to the diffuson and the cooperon. 

8.1 The incoherent albedo: the diffuson 

We have calculated previously the various contributions to the intensity. The incoherent part is 
given by the diffuson contribution, namely the relation (^), where we take as the source term the 
incident plane wave 

^(ri) = y^-l-- \/^e-^ (77) 

where r is the impact point on the surface z = and ri is the position of the first collision. The 
albedo rewrites 

ad = T / dr i dr 2|G B (r2,R)| 2 r(r 1 ,r 2 )e^l ri - r '/'= (78) 

The Fraunhoffer approximation consists in taking the limit |R — r 2 | — > oo, so that the Green 
function G (r 2 , R) can be expanded as 



G (r 2 ,R) = 



-r 2 \/2h 



gifc|R-r 2 | 

4vr|R-r 2 | 



p — ikR 
r 2 |/2Z e -ifcs e .r 2 e 



e-— — (79) 



By defining the respective projections p and p,o of the vectors §i and s e on the z axis, we obtain 
finally 

a d = ,. , ?c , dr 1 dr 2 e "o'^e ^r(r 1} r 2 ) (80) 



This expression calls for a number of remarks. We first notice that it does not contain any 
dependence on the direction namely on the incident and emerging vectors §i and s e . Therefore, 
at this approximation, the albedo is flat i.e. does not have any angular structure. Second, if we 
compare this expression to its counterpart for the conductivity, we see that although both of them 
describe transport, the conductivity vanishes for the diffuson while it does not for the albedo. 
The structure factor T is related to the diffusion probability through the relation 

Pd(r 1 ,r a ) = -/s-r(n,r a ) (81) 
Aire 



so that we obtain 



4tt/2 



a d = ~Tr? I dz 1 dz 2 e e We I d 2 pP d (zi, z 2 , p) (82) 



s 



where now due to the geometry of the medium, the function Pd depends on zi, z 2 and the 
projection p of the vector ri — r 2 onto the plane z = 0. Within the diffusion approximation, Pd 
can be calculated for this geometry using the image method flf @. This gives 

where r'* = (r'*i,—z' + 2z ) and z = |Z e . Then, it appears that Pd does not vanish on the 
plane z = but instead for z = — z . ft is possible for that geometry to compare the validity of 
the diffusion approximation with the known exact solution for the semi-infinite problem without 
sources (the so called Milne problem) that can be solved using the Wiener-Hopf method [ p0| . 
There, we obtain instead that Pd vanishes on the plane z = — 0.7104Z e . This justifies the use of 
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the diffusion approximation. But, it is important at this point to make two important remarks. 
First, it would be a mistake to believe that on the basis of the two examples, infinite and semi- 
infinite spaces, for which the diffusion approximation works well that it is indeed always the case. 
The generalization and the validity of the diffusion equation in a restricted geometry is a difficult 
problem. Second, all the expressions we have obtained so far are valid for isotropic scatterers for 
which there is no difference between the elastic and the transport mean free paths. When relaxing 
this approximation, namely dealing with anisotropic scattering, we have to include both of them. 
Although it appears clearly from simple physical considerations that in the expression (|32|), we 
must leave the elastic mean free path l e in the exponential factors that describe the first and last 
collisions, and replace it by the transport mean free path in the expression of Pa, this point needs 
further investigation using for instance transport theory j(| [l^] . 
Finally, by replacing ( |8^ ) into (|8^), we obtain 

oid = -pMMoi-p + — — ) (84) 

47T l e /! + flo' 

8.2 The coherent albedo: the cooperon 

Along the same lines, we can now evaluate the contribution of the cooperon to the albedo. Using 
the expression (|50| ) for the intensity, we obtain 



a c (s h § e ) = j—^ J dndr 2 e"^- + ™ ^^r(r 1; r 2 )e lfc ^ +s ^-^-r 1 ; (85) 

Unlike the previous case, there is now a phase term present in this expression which is at the origin 
of a new angular structure of the albedo. It is straightforward to check that, assuming §i + s e = 0, 
i.e. measuring the albedo right in the backscattering direction, we have 

a c {9 = 0)=a d (86) 

where the angle 9 is between the directions §i and i e . Then, the total albedo a(8) = ad + a c (9) 
is such that 

a{6 = 0) = 2a d (87) 
Using the expression ( |83| ) for the probability and integrating, we obtain 

a ^°) = -5Z n i ..u.i wi i ...u.i t ( n + 2 7-7— ) 



87T (1 + ^fcjj e )(l + fl k±l e ) fcjJe M + MO ' 

in this expression, = (ki + k e )j^ = fc(§i + s e )^ is the projection onto the plane xOy of the 
vector ki + k e . At small angles, we have k± ~ ^-6, while at large angles, a c (8 — > oo) = 0. Then, 
the coherent contribution to the albedo is finite only in a cone of angular aperture -^r around 
the backscattering direction. By expanding near 9 = we obtain the expression 

a e (fl)^ rf -^ > + Zo) V + 0(fc ± ) 2 (89) 

47T l e 

Thus, the albedo shows a cusp near 9 — namely, its derivative is discontinuous at this point. 



9 Example 3: Spectral correlations 

In the previous two examples, we discussed transport coefficients. They give a description of the 
system when it is connected to another "reference" medium. There is another characterization 
independent of the coupling to another system, which focuses on spectral properties i.e. properties 
of the energy spectrum of the solutions of the wave equations. For electronic systems, they are 
related to the equilibrium thermodynamic properties like the magnetization of a metallic sample. 
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For disordered systems, it also raises an important issue: these systems are part of the larger 
class of complex systems, usually non integrable. The complete set of correlation functions of the 
energy levels provides a complete description of the thermodynamic properties. For instance, it 
is possible just by inspection of the energy spectrum to determine whether or not an electronic 
system is a good metal or an Anderson insulator. 

Let us start by defining the set of eigenenergies e a of the Schrodinger (or Hclmholtz) equation 
in a confined geometry (e.g. a box). The density of states per unit volume is given by 

a 

It can be written as well in terms of the Green function for the Schrodinger equation, using the 
equality 



-— Im / drG R {r,r,e) (91) 
7rS2 J 



By averaging over the disorder, we define the various correlation functions. For instance, the 
two-point correlation function is 

p(ei)p(e 2 )-7J(ei) p(e 2 ) (92) 

where the average density of states coincides with its free value po and therefore is related to the 
mean level spacing A through p = . It is important at this point to make the following remark. 
The average we consider here is over the random potential (^J). We could have taken as well the 
average over different parts of the spectrum as is done for instance in the quantum description of 
classically chaotic billiards where there is no random potential to average over. In the diffusion 
approximation, it has been shown numerically that these two ways to average are equivalent [Q, [2l[] 
but this not need to be the general case. We then define the dimensionless correlation function 



K(u) = P(£) f - " } - 1 (93) 
p(e) 2 

It can be expressed in terms of the Green function under the form 

K(u) = A 2 J drdr'K(r, r', w) (94) 

with 

1 



K (r, r', lu) = — ReG R (r, r, e)G A (r', r', e - u) (95) 

where — c represents the cumulant average. The diffuson and cooperon contributions to the local 
function K(r,r',uj) can be expressed in terms of P d and the structure factor as 

K d (r, r', u) = ^2 Re [P d (v, r'u)P d {r', r, u)] , (96) 

while 

K c (r,r',u;) = 2pl^Re[Ti(v,v')Ti(v',v)} (97) 

when the system has time-reversal invariance. Finally, we define the Fourier transform K(t) which 
is often called the form factor. By collecting the two previous contributions, it can be written 

K{t) = ^\t\Z{\t\) = ^\t\ J p d ( T ,T,\t\)dr (98) 

which is precisely the form that is obtained assuming the Random Matrix Theory |l| provided 
we consider the regime for which h/A ^> t ^>> where td — L 2 /D. The first inequality which 
involves the Hciscnbcrg time h/A enforces the condition of a continuous spectrum while the second 
one indicates that we must be in the ergodic limit where the diffusing particle explored the whole 
system. 
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10 Conclusion 



We have presented in this short review a very partial selection of highlights in the field of multiple 
scattering of waves in disordered media. Our aim was more to give a feeling of the profound unity 
of the physical phenomena and therefore of the methods to handle them rather than discussing an 
extensive list of effects. Among them, we should mention the study of fluctuations of the intensity 
for waves both in the weak disorder regime and near the Anderson transition that will be covered 
in details by the reviews of R. Pnini, A. Genack and G.Maret. This is certainly a problem for 
which the recent developments have been very spectacular. We have studied the coherent albedo 
in the weak disorder case. It should be extended to the strong disorder limit. Finally, our last 
example on spectral correlations should be considered as another way to obtain the semiclassical 
results reviewed in the contribution of D. Delande. 
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